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We present, in detail, a search for the standard model Higgs boson, H , in final states with a 
charged lepton (electron or muon), missing energy, and two or more jets in data corresponding to 
9.7 fb~^ of integrated luminosity collected at a center of mass energy of •Js = 1.96 TeV with the 
DO detector at the Fermilab Tevatron pp Collider. The search uses b-jet identification to categorize 
events for improved signal versus background separation and is sensitive to associated production 
of the H with a W boson, WH Ivbh; gluon fusion with the Higgs decaying to W boson pairs, 
H — >■ WW — Ivjj; and associated production with a vector boson where the Higgs decays to W 
boson pairs, VH — > VWW — >■ Ivjjjj production (where V = W or Z). We observe good agreement 
between data and expected background. We test our method by measuring WZ and ZZ production 
with Z ^ bb and find production rates consistent with the standard model prediction. For a Higgs 
boson mass of 125 GeV, we set a 95% C.L. upper limit on the production of a standard model Higgs 
boson of 5.8xa"sM, where ctsm is the standard model Higgs boson production cross section, while the 
expected limit is 4.7 x asm- We also interpret the data considering models with fourth generation 
fermions, or a fermiophobic Higgs boson. 

PACS numbers: 14.80.Bn, 13.85.Rm 



I. INTRODUCTION 



The Higgs boson is the massive physical state that 
emerges from electroweak symmetry breaking in the 
Higgs mechanism [1-3]. This mechanism generates the 
masses of the weak gauge bosons and explains the fermion 
masses through their Yukawa couplings to the Higgs bo- 
son field. The mass of the Higgs boson {Mh) is a free pa- 
rameter in the standard model (SM). Precision measure- 
ments of various SM electroweak parameters constrain 
Mh to be less than 152 GeV at the 95% C.L. [4-6]. Direct 
searches at the CERN e+e" Collider (LEP) [7] exclude 
Mh < 114.4 GcV at the 95% C.L. The ATLAS and CMS 
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Collaborations, using pp collisions at the CERN LHC, ex- 
clude masses from 110 < Mh < 600 GcV, except for a 
narrow region between 122 and 127 GcV [8, 9]. Both ex- 
periments observe a resonance at a mass of w 125 GeV, 
primarily in the 77 and ZZ final states, with a sig- 
nificance greater than 5 standard deviations (s.d.) that 
is consistent with SM Higgs boson production [10, 11]. 
The CDF and DO Collaborations at the Fermilab Teva- 
tron Collider report a combined analysis that excludes 
the region 147 < Mh < 179 GeV [12] and shows ev- 
idence at the 3 s.d. level for a particle decaying to bb, 
produced in association with a, W or Z boson, consis- 
tent with SM WH/ZH production [13]. Demonstrating 
that the observed resonance is the SM Higgs boson re- 
quires also observing it at the predicted rate in the bh 
final state, which is the dominant decay mode for masses 
below Mh < 135 GeV. 

The dominant production process for the Higgs bo- 
son at the Tevatron Collider is gluon fusion {gg H), 
followed by the associated production of a Higgs boson 
with a vector boson (VH), then via vector boson fusion 
{VVqq' Hqq'). At masses below Mh « 135 GeV, the 
Higgs boson mainly decays to a pair of b quarks, while 
for larger masses, the dominant decay is to a pair of W 
bosons. Because the H ^ bb process is difficult to dis- 
tinguish from background at hadron colliders, it is more 
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effective to search for tlic Higgs boson produced in asso- 
ciation witli a vector boson for this decay channel. 

This Article presents a search by the DO collabora- 
tion for the SM Higgs boson using events containing one 
isolated charged lepton {£ = e or /i), a significant im- 
balance in transverse energy (^t), and two or more jets. 
It includes a detailed description of the WH — >■ li/bb 
search, initially presented in Ref. [14] and used as an in- 
put to the result presented in Ref. [13], differing from 
and superseding that result due to an updated treat- 
ment of some systematic uncertainties as described in 
Sec. X below. The complete analysis comprises searches 
for the production and decay channels: WH — ii'bb, 
H WW* ii^jj (where j = u,d,s,c), and VH 
VWW* — ivjjjj (where V = W or Z). This search also 
considers contributions from ZH production and from 
the decay H ZZ when one of the charged leptons 
from Z —i' ££ decay is not identified in the detector. We 
optimize the analysis by subdividing data into mutually 
exclusive subchannels based on charged lepton flavor, jet 
multiplicity, and the number and quality of candidate 
b quark jets. This search also extends the most recent 
DO WH — >■ £i/bb search [14] by adding subchannels with 
looser 6-quark jet identification requirements and sub- 
channels with four or more jets. These additional sub- 
channels are primarily sensitive to — > WW* — )- £vjj 
and VH — > VWW* — ?> £vijii production and extend 
the reach of our search to Mh = 200 GeV. We present a 
measurement of V Z production with Z — >■ 66 as a cross 
check on our methodology in Sec. XI. In addition to 
our standard model interpretation, we consider interpre- 
tations of our result in models with a fourth generation 
of fcrmions, and models with a fermiophobic Higgs as 
described in Sec. XIII. 

Several other searches for WH — > £vbb production 
have been reported at a pp center-of-mass energy of 
a/s = 1.96 TeV, most recently by the CDF Collabora- 
tion [15]. The results presented here supersede previous 
searches by the DO Collaboration, presented in Refs. [16- 
20], which used subsamples of the data presented in this 
Article. They also supersede a previous search for Higgs 
boson production in the £vii final state by the DO Col- 
laboration [21]. 



II. THE DO DETECTOR 

This analysis relies on all major components of the 
DO detector: tracking detectors, calorimeters, and muon 
identification system. These systems are described in 
detail in Ref. [22-25]. 

Closest to the interaction point is the silicon microstrip 
tracker (SMT) followed by the central scintillating fiber 
tracker (CFT) . These detector subsystems are located in- 
side a 2 T magnetic field provided by a superconducting 
solenoid. They track charged particles and are used to 
reconstruct primary and secondary vertices for pseudo- 
rapidities [26] of |?7| < 3. Outside the solenoid is the liq- 



uid argon/uranium calorimeter consisting of one central 
calorimeter (CC) covering |?7| ^ 1 and two end calorime- 
ters (EC) extending coverage to I77I « 4. Each calorime- 
ter contains an innermost finely segmented electromag- 
netic layer followed by two hadronic layers, with fine and 
coarse segmentation, respectively. The main functions of 
the calorimeters arc to measure energies and help identify 
electrons, photons, and jets using coordinate information 
of significant energy clusters. They also give a measure 
of the l^T- A preshower detector between the solenoidal 
magnet and central calorimeter consists of a cylindrical 
radiator and three layers of scintillator strips covering the 
region \ri\ < 1.3. The outermost system provides muon 
identification. It is divided into a central section that 
covers |ry| < 1 and forward sections that extend coverage 
out to \rj\ « 2. The muon system is composed of three 
layers of drift tubes and scintillation counters, one layer 
before and two layers after a 1.8 T toroidal magnet. 

III. EVENT TRIGGER 

Events in the electron channel are triggered by a logi- 
cal OR of triggers that require an electromagnetic object 
and jets, as described in Ref. [20]. Trigger efficiencies are 
modeled in the Monte Carlo (MC) simulation by apply- 
ing the trigger efficiency, measured in data, as an event 
weight. This efficiency is parametrized as a function of 
electron 77, azimuthal angle (j) [27], and transverse mo- 
mentum pt- For the events selected in our analysis, these 
triggers have an efficiency of (90 — 100)% depending on 
the trigger and the region of the detector. 

The muon channel uses an inclusive trigger approach, 
based on the logical OR of all available triggers, ex- 
cept those containing lifetime-based requirements that 
can bias the performance of 6-jet identification. To de- 
termine the trigger efficiency, we compare data events se- 
lected with a well-modeled logical OR of the single muon 
and muon-|-jets triggers (T^or)j which are about 70% ef- 
ficient, to events selected using all triggers. The increase 
in event yield in the inclusive trigger sample is used to de- 
termine an inclusive trigger correction for the MC trigger 
efficiency, Pcorr, relative to the T^oR trigger ensemble: 

_ (A^Data — Nyi3)incl — (A^Data " -^Mj)t^or /-, \ 

where the numerator is the difference between the num- 
ber of data events in the inclusive trigger sample and the 
Tp,oR trigger sample, after subtracting off instrumental 
multijet (MJ) backgrounds, and the denominator is the 
number of MC events (after the event selection and nor- 
malization to data described in Sec. VIII and the MC 
corrections are applied as described in Sec. VIA) with 
the trigger efficiency set to 1 . The total trigger efficiency 
estimate for events in the muon channel is 7)iOR + ^corr, 
limited to be < 1. 

Triggers based on jets and make the most signifi- 
cant contributions to the inclusive set of triggers beyond 
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those included in the well- modeled T^oR trigger set. To 
account for these contributions, the correction from T^oR 
triggers to the inclusive set of triggers is parametrized as 
a function of the scalar sum of the transverse momenta 
of all jets, Ht^ and the ^t, and is derived for separate 
regions in muon 77. 

For I?/! < 1.0, events are dominantly triggered by sin- 
gle muon triggers, while for \ri\ > 1.6, triggers based 
on the logical OR of muon-|-jets prevail. The third re- 
gion, 1.0 < I77I < 1.6, is a mixture of single muon and 
muon-|-jets triggers. In the \r]\ < 1.0 and 1.0 < jryj < 1.6 
regions the detector support structure allows only partial 
coverage by the muon system. This impacts the muon 
trigger efficiency in the region — 2<0<— 1.2. In these 
regions, we therefore derive separate corrections. The in- 
clusive trigger approach results in a gain of about 30% in 
efficiency over using only muon and the muon-t-jets trig- 
gers. Examples of these corrections, Pcorr are shown in 
Fig. 1. 
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FIG. 1: (color online) Data-derived muon trigger correction 
to account for the resulting efficiency gain in moving from 
single muon and muon-|-jets triggers to inclusive triggers as 
a function of Ht for < 1.0, shown (a) for events with 
< 50 GeV and (b) for events with > 50 GeV. The 
black circles show the correction when the muon is in the 
region of (j) (—2 < (j) < —1.2) where there is a gap in the 
muon coverage for detector supports, and the red triangles 
show the correction elsewhere in S. 



IV. IDENTIFICATION OF LEPTONS, JETS, 
AND 

To reconstruct the candidate W(— >■ iv) boson, our se- 
lected events are required to contain a single identified 
electron or muon together with significant ■ To ensure 
statistical independence with channels that contain more 
than one lepton, we do not consider events with more 
than one electron or muon. Two or more jets are also re- 
quired in order to study WH — )- Ivhh, H — >■ WW — >■ tvjj, 
and VH — >■ VWW — >■ Ivjjjj production. Two sets of 
lepton identification criteria are applied for each lepton 
channel in order to form a "loose" sample, used to es- 
timate the multijet background from data as described 
in Sec. VII, and a "tight" sample used to perform the 
search. The event selection procedure, prior to 6-jet cat- 
egorization, is similar to that described in Ref. [20] and 
described in more detail below. 

Electrons with > 15 GeV are selected in the pseu- 
dorapidity regions \ri\ < 1.1 and 1.5 < \ri\ < 2.5, corre- 
sponding to the CC and EC, respectively. A multivariate 
discriminant is used to identify electrons. The discrimi- 
nant is based on a boosted decision tree [28-32] (BDT) as 
implemented in the tmva package [33] with input vari- 
ables that are listed below. The BDTs are discussed in 
more detail in Sec. IX. The loose and tight electron sam- 
ples are defined by different requirements on the response 
of this multivariate discriminant that are chosen to re- 
tain high electron selection efficiencies while suppressing 
backgrounds at differing rates. 

Leptons coming from the leptonic decays of W bosons 
tend to be isolated from jets. Isolated electromagnetic 
showers are identified within a cone in rj-cf) space of ATZ = 
y/Ar]^ + A02 < 0.4 [34]. In the CC (EC), an electromag- 
netic shower is required to deposit 97% (90%) of its total 
energy within a cone of radius ATZ ~ 0.2 in the electro- 
magnetic calorimeter. The showers must have transverse 
and longitudinal distributions that are consistent with 
those expected from electrons. In the CC region, a re- 
constructed track, isolated from other tracks, is required 
to have a trajectory that extrapolates to the electromag- 
netic (EM) shower. The isolation criteria restrict the sum 
of the scalar pT of tracks with pT > 0.5 GeV within a 
hollow cone of radius 0.05 < ATZ < 0.4 surrounding the 
electron candidate to be less than 2.5 GeV. The BDT 
is constructed using additional information such as: the 
number and scalar px sum of tracks in the cone of radius 
ATZ < 0.4 surrounding the candidate cluster, track-to- 
cluster-matching probability, the ratio of the transverse 
energy of the cluster to the transverse momentum of the 
track associated with the shower, the EM energy frac- 
tion, lateral and longitudinal shower shape characteris- 
tics, as well as the number of hits in the various layers of 
the tracking detector, and information from the central 
preshower detector. The discriminants are trained using 
Z/7* — > ee data events. 

We select muons with pT > 15 GeV and \r]\ < 2.0. 
They are required to have reconstructed track segments 
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in layers of the muon system both before and after the 
toroidal magnet, except where detector support structure 
limits muon system coverage, for which the presence of 
track segments in any layer is sufficient. The local muon 
system track must be spatially matched to a track in the 
central tracker. 

Muons originating from semi-leptonic decays of heavy 
flavored hadrons are typically not isolated due to jet 
fragmentation and secondary particles from the partial 
hadronic decays. We employ a loose muon definition, re- 
quiring minimal separation of ATZ{fi,j) > 0.5 between 
the muon and any jet, while the tight identification has 
additional isolation requirements. For tight muons, the 
scalar sum of the pt of tracks with ATZ < 0.5 around the 
muon candidate is required to be less than 0.4 x p^. Fur- 
thermore, the transverse energy deposits in the calorime- 
ter in a hollow cone of 0.1 < ATZ < 0.4 around the muon 
must be less than 0.12 x pj.. To suppress cosmic ray 
muons. scintillator timing information is used to require 
hits in the detector to coincide with a beam crossing. 

To reduce backgrounds from Zf-f* — ^£-|-jcts and ti 
production, we reject events containing more than one 
tight-isolated charged lepton. 

Jets are reconstructed in the calorimeters in the region 
\r]\ < 2.5 using an iterative midpoint cone algorithm, with 
a cone size of ATZ = 0.5 [35]. To minimize the possibility 
that jets are caused by noise or spurious energy deposits, 
the fraction of the total jet energy contained in the elec- 
tromagnetic layers of the calorimeter is required to be be- 
tween 5% and 95%, and the energy fraction in the coarse 
hadronic layers of the calorimeter is required to be less 
than 40%. To suppress noise, different energy thresholds 
are also applied to clustered and to isolated cells [36]. 
The energy of the jets is scaled by applying a correction 
determined from 7-l-jet events using the same jet-finding 
algorithm. This scale correction accounts for additional 
energy (e.g., residual energy from previous bunch cross- 
ings and energy from multiple pp interactions) that is 
sampled within the finite cone size, the calorimeter en- 
ergy response to particles produced within the jet cone, 
and energy flowing outside the cone or moving into the 
cone via detector effects [36]. We also apply an addi- 
tional correction that accounts for the flavor composition 
of jets [37]. 

Jet energy calibration and resolution are adjusted in 
simulated events to match those measured in data. This 
correction is derived from Z{-^ ee)-|-jet events from the 
Pt imbalance between the Z boson and the recoiling jet in 
MC simulation when compared to that observed in data, 
and applied to jet samples in MC events. Differences in 
reconstruction thresholds in simulation and data are also 
taken into account, and the jet identification efSciency 
and jet resolution are adjusted in the simulation to match 
those measured in data. All selected jets are required to 
have Pt > 20 GeV and [tjI < 2.5. We require that jets 
originate from the primary vertex (PV), such that each 
selected jet is matched to at least two tracks with pT > 
0.5 GeV that have at least one hit in the SMT detector 



and a distance of closest approach with respect to the PV 
of less than 0.5 cm in the transverse plane and less than 
1 cm along the beam axis (2). Interaction vertices are 
reconstructed from tracks that have pt > 0.5 GeV with 
at least two hits in the SMT. The primary vertex is the 
reconstructed vertex with the highest average pT of its 
tracks. Vertex reconstruction is described in more detail 
in Rcf. [38] . We also require that the PV be reconstructed 
within zpv = ±60 cm of the center of the detector. 

The ^T is calculated from individual calorimeter cell 
energies in the electromagnetic and fine hadronic sections 
of the calorimeter and is required to satisfy ^t > 15 GcV 
for the electron channel and ^t > 20 GeV for the muon 
channel. Energy from the coarse hadronic layers that is 
contained within a jet is also included in the ^t calcu- 
lation. A correction for the presence of any muons and 
all energy corrections applied to electrons and jets are 
propagated to the value of ^t ■ 



V. TAGGING OF b-QUARK JETS 

The 6-tagging algorithm for identifying jets originating 
from b quarks is based on a multivariate discriminant us- 
ing a combination of variables sensitive to the presence of 
tracks or secondary vertices displaced significantly in the 
x-y plane from the pp interaction vertex. This algorithm 
provides improved performance over the neural network 
algorithm described in Ref. [38]. 

Jets considered by the 6-tagging algorithm are required 
to be "taggable," i.e., contain at least two tracks with 
each having at least one hit in the SMT. The efficiency 
of this requirement accounts for variations in detector 
acceptance and track reconstruction efficiencies at differ- 
ent locations of the PV prior to the application of the 
6-tagging algorithm, and depends on the z position of 
the PV and the pt and 77 of the jet. For jets that pass 
through the geometrical acceptance of the tracking sys- 
tem, this efficiency is typically about 97%. The efficiency 
for 6-tagging is determined with respect to taggable jets. 
The correction for taggability is measured in the selected 
data sample, while the corrections for 6-tagging are deter- 
mined in an independent heavy-flavor jet enriched sample 
of events that include a jet containing a muon, as de- 
scribed in Ref. [38] . The efficiency for jets to be taggable 
and to satisfy 6-tagging requirements in the simulation is 
corrected to reproduce the respective efficiencies in data. 

We define six independent tagging samples with zero, 
one loose, one tight, two loose, two medium, or two tight 
6-tagged jets. An inclusive "pretag" sample is also con- 
sidered for parts of this analysis. Events with no jets 
satisfying the 6-tagging criteria are included in the zero 
6-tag sample. If exactly one jet is 6-tagged, and the b- 
identification discriminant output for that jet, satis- 
fies the tight selection threshold (fejj-, > 0.15), that event 
is considered part of the one tight 6-tag sample. Events 
with exactly one 6-tagged jet that fails the tight selec- 
tion threshold, but passes the loose selection threshold 
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> 0.02) are included in the one loose 6-tag sample. 
Events with two or more 6-tagged jets are assigned to ei- 
ther the two loose 6-tags, two medium 6-tags, or two tight 
6-tags category, depending on the value of the average b- 
identification discriminant of the two jets with the high- 
est discriminant values, i.e., the double tight category 
is required to satisfy + &jd)/2 > 0.55; the medium 
category is 0.35 < (fej^ + bl^)/2 < 0.55; and the loose 
category is 0.02 < (6f^ 4- b{l)/2 < 0.35 (see Fig. 2). The 
operating point for the loose (medium, tight) threshold 
has an identification efficiency of 79% (57%, 47%) for 
individual b jets, averaged over selected jet px and 77 dis- 
tributions, with a 6-tagging misidentification rate of 11% 
(0.6%, 0.15%) for hght quark jets {If), calculated by the 
method described in Ref. [38]. 
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FIG. 2: (color online) Average of the b-identification discrim- 
inant outputs of each jet in events with two jets. 



VI. MONTE CARLO SIMULATION 

We account for all Higgs boson production and decay 
processes that can lead to a final state containing exactly 
one charged well isolated lepton, ^t, and two or more 
jets. The signal processes considered are: 

• Associated production of a Higgs boson with a vec- 
tor boson where the Higgs boson decays to bb, cc, 
TT, or VV . The associated weak vector boson 
decays leptonically in the case of — >■ 66 and 
either leptonically or hadronically in the case of 
H WW. Contributions from Z{-^ ei)H{~^ bb) 
production arise from identifying only one charged 
lepton in the detector, with the other contributing 
to the ^T- 

• Higgs boson production via gluon fusion with the 
subsequent decay H — > VV, where one weak vector 
boson decays leptonically (with exactly one identi- 
fied lepton). 



• Higgs boson production via vector boson fusion 
with the subsequent decay H — >■ VV, where one 
weak vector boson decays leptonically (with exactly 
one identified lepton). 

Various SM processes can mimic expected signal sig- 
natures, including F-|-jets, diboson (VV), MJ, tt, and 
single top quark production. 

All signal processes and most of the background pro- 
cesses are estimated from MC simulation, while the 
MJ background is evaluated from data, as described in 
Sec. VH. We use pythia [39] to simulate all signal pro- 
cesses and diboson processes. The X^-f-jets and samples 
are simulated with the ALPGEN [40] MC generator in- 
terfaced to PYTHIA for parton showering and hadroniza- 
tion, while the SINGLETOP event generator [41, 42] inter- 
faced to PYTHIA is used for single top quark events. To 
avoid overestimating the probability of further partonic 
emissions in pythia, the MLM factorization ("match- 
ing") scheme [43] is used. All of these simulations use 
CTEQ6L1 [44, 45] parton distribution functions (PDF). 

A full GEANT-based [46] detector simulation is used 
to process signal and background events. To account 
for residual activity from previous beam crossings and 
contributions from the presence of additional pp interac- 
tions, events from randomly selected beam crossings with 
the same instantaneous luminosity profile as the data are 
overlaid on the simulated events. All events are then re- 
constructed using the same software as used for data. 

The signal cross sections and branching fractions are 
normalized to the SM predictions [12]. The WH and ZH 
cross sections are calculated at next-to-next-to-leading 
order (NNLO) [47], with MSTW2008 NNLO PDFs [48]. 
The gluon fusion process uses the NNLO-|-next-to-next- 
to- leading log (NNLL) calculation [49], and the vector 
boson fusion process is calculated at NNLO in QCD [50]. 
The Higgs boson decay branching fractions are obtained 
with HDECAY [51, 52]. We use NLO cross sections to 
normalize single top quark [53] and diboson [54, 55] pro- 
duction, while we use an approximate NNLO calculation 
for tt production [56] . The pr of the Z boson in Z+jcts 
events is corrected to match that observed in data [57]. 
The pt of the W boson in VF-|-jets events is corrected us- 
ing the same dependence but taking into account the dif- 
ferences between the pr spectra of the Z and W bosons in 
NNLO QCD [58]. Additional scale factors to account for 
higher order terms in the ALPGEN MC for the T^-l-heavy 
fiavor jets, V + hf, are obtained from mcfm [55, 59]. The 
T^-|-jcts processes are then normalized to data for each 
lepton fiavor and jet multiplicity separately as described 
in Sec. VIII. 



A. MC Reweighting 

Motivated by previous comparisons of ALPGEN with 
data [60] and with other event generators [43], we de- 
velop corrections to VF -|- jets and Z + jets MC samples 
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to correct for the shape discrepancies in kinematic dis- 
tributions between data and simulation. The corrections 
are derived based on the direct comparison between data 
and MC samples prior to the application of 6-tagging, 
where any contamination from signal is very small. 

To improve the description of jet directions, we correct 
the r] distributions of the leading and second leading jets 
in W/Z+}cts events. The correction function is a fourth- 
order polynomial determined from the ratio of the l^-f-jets 
events in MC and data minus non-V^-l-jets backgrounds. 
The modeling of the lepton rj in PF-|-jets events is ad- 
justed by applying a second-order polynomial correction. 
Correlated discrepancies observed in the leptonically de- 
caying W boson transverse momentum, , and the jet 
angular separation, AT?, (j 1,^2), arc corrected through 
two reweighting functions in the two-dimensional ATZ- 
plane [20]. The reweighting is applied only to 
W + jets events, while the ATZ reweighting is applied to 
both -I- jets and Z -|- jets events. Each of these correc- 
tions is designed to change differential distributions, but 
to preserve normalization. Corrections are on the order 
of a few percent in the highly populated region of each 
distribution and may exceed 10% for extreme values of 
each distribution. 

All corrections are derived in events selected with 
muon-fjets triggers to minimize uncertainties due to con- 
tamination from MJ events, and are applied to both the 
electron and muon channels. Additional Pj^, ATZ, and 
lepton 77 corrections and corresponding systematic uncer- 
tainties are determined from events selected with inclu- 
sive muon triggers and are applied to events containing 
muons, accounting for variations in modeling distribu- 
tions of the inclusively triggered events. 



VII. MULTIJET BACKGROUND 

The MJ background, events where a jet is misidcntified 
as a lepton, is determined from the data prior to the ap- 
plication of 6-tagging, using a method similar to the one 
used in Ref. [20]. This method involves applying event 
weights that depend on the relative efBciency ef^p, of a 
lepton passing loose requirements to subsequently pass 
the tight requirements and a similar relative probability, 
P]^'', for a MJ event to pass these sequential selections. 
A MJ template is constructed by selecting events from 
data in which the lepton passes the loose isolation re- 
quirement, but fails the tight requirement, as described 
in Sec. IV. Each event in the MJ template is weighted 

by 



pMJ 
1 _ pMJ ' 



(2) 



where Pj^"' is a function of the event kinematics. Since 
the MJ template contains a contribution from events 
with real leptons originating from Icptonic decays of W/Z 
bosons, we correct the normalization of the V-f jets MC 



using the event weight 



pMJ 



-LTJ 



(3) 



where efx and are functions of event kinematics. 

The efficiencies efj^ are functions of lepton px, and they 
are determined from Z/j* —>■ U events. The probabili- 
ties are determined in the region 5 < < 15 GeV 
from the measured ratio of the number of events with 
tight leptons and those with loose leptons after correct- 
ing each sample for the expected MC contribution from 
real leptons in the specific kinematic interval. Electron 
channel probabilities are parametrized in pT , calorimeter 
detector 77, and rnXn A(p{]^T , j) while probabilities in the 
muon channel are parametrized in pT for different regions 
in muon detector 77 and A(/)(^t,/^)- 



VIII. EVENT SELECTION 



Events are required to have one isolated charged lep- 
ton, large and two or more jets, as described in 
Sec. IV. To suppress MJ backgrounds, events must sat- 
isfy the additional requirement that Mjf > 40 GeV — 
0.5 X ]^T^ where is the transverse mass [61] of the 
W boson. We then perform the final normalization of 
the T^-|-jets and MJ backgrounds via a simultaneous fit to 
data in the MjT distribution after subtracting the other 
SM background predictions from the data as described in 
Refs. [14, 19, 20]. The distribution of M]f after this nor- 
malization procedure is shown in Fig. 3(a). We perform 
separate fits for each lepton flavor and jet multiplicity 
category before dividing events into categories based on 
the number and quality of identified b jets, as described 
in Sec. V. All events passing these selection criteria con- 
stitute the pretag sample, and each pretag event also 
belongs to exactly one of the six independent 6-tag cate- 
gories. Only the zero and one-loose 6-tag categories are 
considered when searching for the signal in events with 
four or more jets because ti production dominates the 
small amount of signal present in higher 6-tag categories. 

The expected number of events from each signal and 
background category is compared to the observed data 
for each 6-jet identification category for events with two 
jets, three jets, and four or more jets in Tables I, II, 
and III, respectively. Selected kinematic distributions are 
shown for all selected events in Figs. 3 and 4, and the di- 
jet invariant mass for events with two jets is shown for all 
&-tag categories in Figs. 5 and 6. In all plots, data points 
are shown with error bars that rcfiect the statistical un- 
certainty only. Discrepancies in data-MC agreement are 
within our systematic uncertainties described in Sec. X. 
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TABLE I: Observed number of events in data and expected number of events from each signal and background source (where 
V — W, Z) for events with exactly two jets. The expected signal is quoted at Mr = 125 GeV. The total background uncertainty 
includes all sources of systematic uncertainty added in quadrature. 





Pretag 


6-tags 


1 loose b-tag 


1 tight b-tag 


2 loose b-tags 


2 med. b-tags 2 tif 


;ht b-tags 


VH -5> lubb 


37.3 


6.4 


4.0 


11.6 


3.2 


4.6 


7.7 


H ^VV ^ lujj 


24.7 


18.8 


3.9 


1.8 


0.3 


0.07 





VH VVV -)■ tvjjjj 


13.0 


9.3 


2.3 


1.2 


0.3 


0.04 


0.01 




ouou 




yuo 




1 no 




9S 
oo 


V + (p,u,d, s)-jets 


182 271 


148 686 


26 421 


6174 


1762 


132 


13 


V + {hi Ice) 


27443 


15 089 


4872 


5236 


978 


691 


691 


top {tt + single top) 


3528 


758 


455 


1289 


247 


333 


462 


Multijet 


58 002 


43 546 


9316 


3700 


946 


298 


195 


Total expectation 


276 930 


212114 


42 032 


16 935 


4043 


1496 


1400 


Total uncertainty 


± 14 998 


± 11352 


± 2438 


± 1696 


± 362 


± 117 


± 175 


Observed events 


276 929 


211169 


42 774 


16 406 


4057 


1358 


1165 



TABLE 11: Observed number of events in data and expected number of events from each signal and background source (where 
V = W,Z) for events with exactly three jets. The expected signal is quoted at AIh ~ 125 GeV. The total background 
uncertainty includes all sources of systematic uncertainty added in quadrature. 





Pretag 


b-tags 


1 loose b-tag 1 tif 


;ht b-tag 


2 loose b-tags 2 med. 


b-tags 2 tif 


;ht b-tags 


VH ivbh 


8.6 


1.3 


1.0 


2.4 


0.9 


1.1 


1.7 


H ^VV ^ £ujj 


8.8 


6.0 


1.7 


0.8 


0.3 


0.07 


0.01 


VH VVV -> iujjjj 


7.3 


4.5 


1.6 


0.9 


0.3 


0.05 


0.01 


Diboson 


1138 


727 


238 


113 


42 


14 


10 


V + (fir, u, d, s)-jets 


24 086 


18 078 


4577 


976 


582 


34 


3 


V + {bb/cc) 


6625 


3213 


1349 


1250 


411 


228 


164 


top {tt + single top) 


3695 


563 


419 


1123 


365 


460 


570 


Multijet 


10 364 


6629 


2162 


933 


367 


130 


82 


Total expectation 


45 908 


29 209 


8746 


4395 


1768 


867 


830 


Total uncertainty 


± 2582 


± 1619 


± 587 


± 528 


± 209 


± 118 


± 113 


Observed events 


45 907 


28 924 


8814 


4278 


1815 


879 


797 



TABLE III: Observed number of events in data and expected 
number of events from each signal and background source 
(where V = W, Z) for events with four or more jets. The 
expected signal is quoted at Mh = 125 GeV. The total back- 
ground uncertainty includes all sources of systematic uncer- 
tainty added in quadrature. 





Pretag b-tags 


1 loose b-tag 


VH eubb 


1.4 


0.2 


0.2 


H ^VV ^ lujj 


2.4 


1.4 


0.6 


VH ~> VVV -> tvjjjj 


3.6 


2.0 


0.8 


Diboson 


199 


112 


46 


V + {g,u,d,s)-]ets. 


3055 


2143 


679 


V + {bb/cc) 


1280 


542 


286 


top {tt + single top) 


2889 


311 


268 


Multijet 


2092 


1110 


450 


Total expectation 


9516 


4217 


1729 


Total uncertainty 


± 530 


± 231 


± 144 


Observed events 


9685 


3915 


1786 



IX. MULTIVARIATE SIGNAL 
DISCRIMINANTS 

We employ multivariate analysis (MVA) techniques to 
separate signal from background events. To separate sig- 
nal from the MJ events, we use a boosted decision tree 
implemented with the tmva package [33] . This multi- 
variate analysis is described in Sec. IX A. A BDT is also 
used to separate signal from other specific background 
sources in events with four or more jets (see Sec. IX D). 
For the final multivariate analysis, we use a BDT in the 
one tight fe-tag channel and all three two &-tag channels, 
and we use a random forest decision tree (RF) [62] imple- 
mented in the statpatternrecognition (SPR) pack- 
age [28, 63] for events in the zero and one loose 6-tag 
channels. 

The BDT and the RF are forms of machine learn- 
ing techniques known as decision trees. Decision trees 
operate on a series of yes/no splits on events that are 
known to be classified as either signal or background. 
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FIG. 3: (color online) Distributions for all selected events 
with two jets of (a) transverse mass of the lepton-_^T system, 
(b) charged lepton pT, and (c) |?t. The signal is multiplied 
by 1000. Overflow events are added to the last bin. 



FIG. 4: (color online) Distributions for all selected events 
with two jets of (a) leading jet pr, (b) second-leading jet pr, 
and (c) AT?, between the leading and second-leading jets. The 
signal is multiplied by 1000. Overflow events are added to the 
last bin. 



11 



V(^lv)+2 jets, zero b-tags 



35000 
30000 
25000 
20000 
15000 
10000 
5000 



V(^lv)+2 jets, two loose b-tags 



D0, 9.7 fb- 
(a) 



+ Data 

■ v+hf 
V+lf 
B Multijet 
□ Signal (xlOOO) 



50 100150200250300350400 

Dijet mass [GeV] 



D0,9.7 fb"" 




50 100150200250300350400 

Dijet mass [GeV] 



V(^lv)+2 jets, one loose b-tag 



> 12000 
O 

O 10000 
o 

^ 8000 

^ 6000 
0) 

> 4000 
HI 

2000 



D0,9.7 fb- 
(b) 



■ 



+ Data 

B vv 

^Top 
■ v+hf 

V+lf 
B Multijet 
□ Signal (x 500) 

Mh=125 GeV 



50 100150200250300350400 

Dijet mass [GeV] 



> 

0) 

O 
o 

CM 
W 

c 
> 

UJ 



lv)+2 jets, two medium b-tags 




50 100150200250300350400 

Dijet mass [GeV] 



5000 

> 

Q 4000 



V(^lv)+2 jets, one tight b-tag 
D0, 9.7 fb-^ +Data 

(C) "^^ 



3000 



C 2000 

> 

HI 

1000 



I Top 

■ v+hf 
V+lf 

■ lUlultijet 

□ Signal (x 200) 




50 100150200250300350400 

Dijet mass [GeV] 



> 

O 
o 

CM 
W 

c 

> 

UJ 



V( 
400 

350 

300 

250 

200 

150 

100 

50 





^lv)+2 jets, two tight b-tags 
D0, 9.7fb"' +Data 



(C) 



L Top 

■ v+hf 
V+lf 

■ lUlultijet 

I I Signal (x 50) 



50 100150200250300350400 

Dijet mass [GeV] 



FIG. 5: (color online) Invariant mass of the leading and 
second-leading jets in events with two jets and (a) zero b- 
tags, (b) one loose b-tag, and (c) one tight fe-tag. The signal 
is multiplied by 1000, 500, and 200, respectively. Overflow 
events are added to the last bin. 



FIG. 6: (color online) Invariant mass of the leading and 
second-leading jets in events with two jets and (a) two loose 
6-tags, (b) two medium 6-tags, and (c) two tight 6-tags. The 
signal is multiplied by 200, 50, and 50, respectively. Overflow 
events are added to the last bin. 
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The splitting is done to maximally separate signal from 
background. The resulting nodes are continually split to 
optimally separate signal from background until either a 
minimum number of events in a node is reached or the 
events in a node are pure signal or pure background. The 
technique of boosting in the BDT builds up a series of 
trees where each tree is retrained, boosting the weights 
for events that are misclassified in the previous training. 
The RF technique creates a collection of decision trees 
where each tree is trained on a subset of the training 
data that is randomly sampled. 

We train separate BDTs and RFs for each lepton fla- 
vor, jet multiplicity, and tagging category, and for each 
hypothesized Higgs boson mass in steps of 5 GeV. Since 
the branching fraction for the Higgs decay to b quarks 
is only significant over the mass range 90-150 GeV, we 
restrict the search in the one tight and two &-tag chan- 
nels to this range of Mh- In the zero and one loose 6-tag 
channels, the primary signal contribution is from Higgs 
decays to vector bosons, the search is performed over the 
mass range of 100-200 GeV. 

Each of the final BDTs and RFs are trained to dis- 
tinguish the signal from all of the backgrounds. We 
choose variables to train the BDTs and RFs that have 
good agreement between data and background simula- 
tion (since the expected contribution from signal events 
is small) , and so that there is a good separation between 
signal and at least one background. Background and sig- 
nal samples are each split into three independent samples 
for use in training, testing, and performing the final sta- 
tistical analysis with each multivariate discriminant. We 
ensure that the discriminant is not biased towards sta- 
tistical fiuctuations in the training sample by comparing 
the training output to the testing sample. The indepen- 
dent sample used for the limit setting procedure ensures 
that any optimizations performed based on the output 
of the training and testing samples do not bias the final 
limits. 



A. Multivariate multijet discriminators 

We train two separate BDTs to separate the MJ back- 
ground from signal events: one for VH{—> bb, cc, tt) 
signals, MYAujiVH), and one for H VV signals, 
MVAmj {H — >■ VV) . The variables used in training these 
BDTs are chosen to exploit kinematic differences between 
the MJ and signal events, and are documented in Ap- 
pendix A. To improve the training statistics, we com- 
bine signal events for AIh = 120, 125, and 130 GcV 
in training. We find that a BDT trained on this com- 
bination of Higgs boson masses has a similar perfor- 
mance when applied to other masses, eliminating the 
need for a mass dependent MJ discriminant. The BDT 
outputs MYAmj{VH) and MVAmj (i? ^ VV) are shown 
in Fig. 7. The MVAujiVH) and MVAmj (i? ^ VV) dis- 
criminant outputs are used as input variables to the final 
MVAs, as detailed in Appendix A. 



B. Final WH lubb MVA analysis 

In events with two or three jets and one tight 6-tag or 
two 6-tags, the WH — > ivbb process provides the domi- 
nant signal contribution. To separate signal from back- 
ground, we train a BDT on the WH — > Ivbb signal and 
all backgrounds. The lists of input variables to the MVA 
and their descriptions arc included in Appendix A. Fig- 
ures 8 and 9 show examples of some of the most effective 
discriminating variables used in our BDTs for the two- 
jet and three-jet channels, respectively, in the one tight 
6-tag and all two 6-tags channels. Figures 10 and 11 show 
the BDT output for the two and three-jet channels, re- 
spectively, in the one tight 6-tag and all the two 6-tag 
channels. 



C. Final H WW Ivjj MVA analysis 

The H — > WW — s- Ivjj process provides the dominant 
signal in events with two or three jets and zero 6-tags or 
one loose 6-tag, since the W boson decays producing a 
b quark are rare. For signal searches in these channels, 
we apply a multivariate technique based on the RF dis- 
criminant. Events in the above tagging categories are 
examined for 100 < Mh < 150 GcV. Since we do not 
perform the search in the one tight and two 6-tag chan- 
nels for Mh > 150 GeV, events having exactly two or 
three jets in all 6-tagging categories (i.e. pretag events) 
are used in the search for 155 < Mh < 200 GeV. 

To suppress MJ background in the electron channel in 
these subchannels, we select events with MVAmj (^^ ^• 
VV) > -0.4 for Mh < 150 GeV in events with zero 
or one loose tag, and MVAmj{VH) > -0.5 for Mh > 
155 GeV in all events. These requirements were opti- 
mized to maximize the ratio of number of signal events 
to the square root of the number of background events. 
The MJ component in the zero or one loose 6-tag muon 
channel is small, so there is no cut applied to the MJ 
MVA outputs. 

We train a RF on the total signal and background from 
all considered physics processes. We optimize the RF in- 
dependently in the electron and muon channels for each 
6-tag and jet multiplicity category. As the signal shape 
is strongly driven by the signal mass hypothesis, we opti- 
mize the MVA variable list at two different mass points: 
at Mh = 125 GeV for masses below 150 GeV and at 
Mh = 165 GeV for masses above 150 GeV. Because the 
resolution of the reconstructed Higgs boson mass is about 
20 GcV for channels presented in this Article, optimiz- 
ing the input variable list at only these mass points is 
sufficient. Each RF is trained using between 14 and 30 
well modeled discriminating variables formed from kine- 
matic properties of either elementary objects like jets or 
leptons, or composite objects, such as reconstructed W 
boson candidates (see Figs. 12 and 13). The lists of in- 
put variables and their descriptions are included in Ap- 
pendix A. The final RF discriminants for the electron 
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FIG. 7: (color online) The multivariate discriminant output for (a) MVAm3{VH) and (b) MVAmj(JY 
The signal is multiplied by 1000. 



VV), for all events. 



and muon channels are shown in Figs. 14 and 15. 



D. Final VH VWW Ivjjjj MVA analysis 

The majority of signal events with four or more jets 
and zero 6-tags or one loose 6-tag are from the VH — > 
VWW — tvjjii process, but there are significant con- 
tributions from direct production via gluon fusion and 
vector boson fusion. Identification of the Higgs boson 
decay products in VH VWW events is complicated 
by the combinatorics of pairing four jets into two hadron- 
ically decaying vector boson candidates and then two of 
the three total vector boson candidates into the Higgs 
boson candidate. The discriminating variables are dif- 
ferent for fully hadronic and semileptonic Higgs boson 
decays, and determining the Higgs boson candidate for 
an event also determines which of these two decay scenar- 
ios is considered. Variables unique to a particular decay 
scenario are set to a default value outside of the physical 
range of that variable in events reconstructed under the 
alternate decay scenario. To reconstruct the two hadron- 
ically decaying vector boson candidates, we examine the 
leading four jets in an event and choose the jet pairings 
that minimize: 

Eab,cd ^ \mab - Mw\ + \mcd ~ Mw\, (4) 

where mab {mcd) is the invariant mass of the a*^ and b*"^ 
(c*h and d^^) jets, and Mw = 80.4 GeV [71]. The Higgs 
boson candidate is then determined by considering the 
semileptonically decaying W boson and the two hadroni- 
cally decaying vector bosons and selecting the vector bo- 
son candidate pair with the minimum ATZ separation in 
an event, out of the three possible pairings. 

Diverse signal processes contribute to the inclusive 
four-jet channel with relative contributions varying with 



Mh- To help mitigate the effect of having many signal 
and background contributions to this search channel, we 
use two layers of multivariate discriminants to improve 
the separation of signal from background. The first layer 
of training focuses on separating the sum of all signal pro- 
cesses from specific sets of backgrounds. Input variables 
for each background-specific discriminant are selected 
based on the separation power between the total sig- 
nal and the backgrounds being considered. Background- 
specific discriminants are trained to separate the sum of 
all Higgs boson signal processes from three specific back- 
ground categories: ti and single top quark production, 
y-|-jets production, and diboson production. The input 
variables and their descriptions are listed in Appendix A. 
Separate background-specific discriminants are trained 
for each Higgs boson mass point considered. Sample in- 
puts and output responses of the background-specific dis- 
criminants are shown in Figs. 16 and 17, respectively, for 
Mh = 125 GcV. 

The background-specific discriminants are used as in- 
puts to the final RF discriminant that is trained to dis- 
criminate all signal processes from the total background 
contributions. Additional input variables for the final 
discriminant are selected based on their separation power 
between the total signal and the total background, and 
are required to be well modeled. The input variables for 
each lepton and 6-tag category are listed in Appendix A. 
Sample inputs and output responses of the final discrim- 
inants are shown in Figs. 18 and 19, respectively, for 
Mh = 125 GeV. 
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FIG. 8: (color online) Distributions of some of the most significant inputs to the final discriminant in events with exactly two 
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FIG. 14: (color online) Distributions of the final discrimi- 
nant output, after the maximum likelihood fit (described in 
Sec. XII), for events in the following channels: (a) two jets, 
zero 6-tags for Mr ~ 125 GeV, (b) two jets, one loose fe-tag 
for Mh = 125 GeV, (c) two jets, all tags Mh = 165 GeV. 
The signal is multiplied by 500, 500, and 200, respectively. 



FIG. 15: (color online) Distributions of the final discrimi- 
nant output, after the maximum likelihood fit (described in 
Sec. XII), for events in the following channels: (a) three jets, 
zero 6-tags for Mh = 125 GeV, (b) three jets, one loose 6-tag 
for Mh = 125 GeV, (c) three jets, aU tags Mh = 165 GeV. 
The signal is multiplied by 500, 500, and 200, respectively. 
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FIG. 16: (color online) Distributions ol the most significant 
inputs to background-specific multivariate discriminants for 
the > 4-jet subchannels: (a) cos 6{£)e,jCM [69, 72], input to dis- 
criminant against V^-|-jets backgrounds, shown for events with 
zero 6-tags; (b) 5I5jets(^) [73], input to discriminant against 
diboson backgrounds, shown for events with zero 6-tags; (c) 
'Y^PT{t-,^T , ji, j2-, j-i, ji), input to discriminant against top 
quark backgrounds, shown for events with one loose 6-tag. 
The Mh = 125 GeV signal is muhiplied by 250 in (c) and by 
500 in (a) and (b). Overflow events are added to the last bin. 
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FIG. 17: (color online) Distributions of the output of 
background-specific multivariate discriminants for the > 4-jet 
subchannels: (a) discriminant against l/-|-jets backgrounds, 
shown for events with zero 6-tags; (b) discriminant against 
diboson backgrounds, shown for events with zero &-tags; (c) 
discriminant against top quark backgrounds, shown for events 
with one loose fo-tag. The AIh = 125 GeV signal is multiplied 
by 250 in (c) and by 500 in (a) and (b). 



21 



> 
0) 

O 
m 

CM 
(A 

*^ 
c 

> 

UJ 



800 
700 
600 
500 
400 
300 
200 
100 




V(^lv)+4 jets, zero b-tags 



D0, 9.7 fb- 
(a) 



+ Data 

■ v+hf 
V+lf 
B Multijet 
□ Signal (x 500) 




100 200 300 400 500 

rrij , [GeV] 

'l2 



V(^lv)+4 jets, one loose b-tag 



in 
d 

w 
c 

> 

UJ 



500 
400 
300 
200 
100 



D0, 9.7 fb- 
(b) 



+ Data 

Bvv 

^Top 
■ v+hf 

V+lf 
B Multijet 
□ Signal (x 250) 

MH=125GeV 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

A RCljg) 



d 


1800 r 






w 


1600 r 


■*— ' 
c 


1400 r 


> 




LLl 


1200 r 




1000 r 




800 - 




600 r 




400 _ 




200 r 











CO 

d 

v> 
■*— ' 

c 
> 

UJ 



V(^lv)+4 jets, zero b-tags 

D0, 9.7 fb-^ + Data 



(a) 



■ vv 

■ Top 

■ v+hf 
V+lf 

■ Multijet 

□ Signal (x500) 




0.4 0.6 0.8 1 

Final Discriminant 
V(^lv)+4 jets, one loose b-tag 

D0, 9.7 fb"^ + Data 



1000 




0.4 0.6 0.8 

Final Discriminant 



FIG. 18: (color online) Distributions of the most significant 
inputs, other than background-specific multivariate discrimi- 
nants, to the final multivariate discriminants for the > 4-jet 
subchannels: (a) nij^r^i, shown for events with zero fe-tags; 
(b) A7?.(^,j3), shown for events with one loose b-tag. The 
Mh = 125 GeV signal is multiplied by 500 in (a) and 250 in 
(b). Overflow events are added to the last bin. 



FIG. 19: (color online) Distributions of the final discrim- 
inant output, after the maximum likelihood fit (described 
in Sec. XII), at Mh = 125 GeV for the four or more jets 
channels with: (a) zero 6-tags and (b) one loose 6-tag. The 
Mh = 125 GeV signal is multiphed by 500. 
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X. SYSTEMATIC UNCERTAINTIES 

We assess systematic uncertainties on signals and back- 
grounds for each of the jet multiphcity and 6-tag channels 
by repeating the full analysis after varying each source of 
uncertainty by ±1 s.d. We consider uncertainties that af- 
fect both the normalizations and the shapes of our MVA 
outputs. 

We include theoretical uncertainties on the tt and sin- 
gle top quark production cross sections (7% [53, 56]), 
diboson production cross section (6% [54]), V + If pro- 
duction (6%), and V + hf production (20%, estimated 
from MCFM [55, 59]). Since the l/-fjcts experimental 
scaling factors for the three- and four-jet channels are 
different from unity, we apply an additional systematic 
uncertainty on the V^-|-jets samples that is uncorrelated 
across jet multiplicity and lepton flavor bins. The size 
of this uncertainty is taken as the uncertainty from the 
F-l-jets fit to data, described in Sec. VII. 

An uncertainty on the integrated luminosity 
(6.1% [74]) affects the normalization of the expected 
signal and simulated backgrounds. Uncertainties that 
affect the final MVA distribution shapes include jet tag- 
gability (3% per jet), 6-tagging efficiency (2.5%-3% per 
heavy-quark jet), the light-quark jet misidentification 
rate (10% per jet), jet identification efficiency (5%), and 
jet energy calibration and resolution (varying between 
5% and 15%, depending on the process and channel) 
as described in Ref. [20]. We also include uncertain- 
ties from modeling that affect both the shapes and 
normalizations of the final MVA distributions. These 
include an uncertainty on the trigger efficiency in the 
muon channel as derived from the data (3%-5%), lepton 
identification and reconstruction efficiency (5%-6%), the 
MLM matching [40] applied to X^+light-flavor events 
(w 0.5%), the ALPGEN rcnormalization and factorization 
scales, and the choice of parton distribution functions 
(2%) as described in Ref. [20]. The trigger uncertainty 
in the muon channel is calculated as the difference 
between applying a trigger correction calculated using 
the ALPGEN rewcightings derived on the T^oR trigger 
sample and applying the nominal trigger correction. 
Since we reweight our ALPGEN samples, we include 
separate uncertainties on each of the five functions used 
to apply the reweighting. The adjusted functions are 
calculated by shifting the parameter responsible for 
the largest shape variation of the fit by ±1 s.d. then 
calculating the remaining parameters for the function 
using the covariancc matrix obtained from the functional 
fit. 

We determine the uncertainty on the MJ background 
shape by relaxing the requirement from Sec. VIII on Mjf 
to Af^ > 30 GeV — 0.5 x and repeating the analysis 
with this selection in place. The positive and negative 
variations are taken to be symmetric. The uncertainty 
on the MJ rate is 15% (20%) for the electron (muon) 
channel. Since our MJ sample is statistically limited, we 
do not correlate the uncertainties on the rate and shape 



across the subchannels. Since we simultaneously fit MJ 
and y-|-jets to match data, we apply a normalization 
uncertainty to the T^-l-jets samples that is anticorrelated 
with the MJ normalization systematics and scales as the 
relative MJ to y-|-jets normalization. 



XI. WZ AND ZZ PRODUCTION WITH Z ^ bb 

_The SM processes W{-> e.v)Z{-^ hh) and Z{-^ U)Z{-^ 
bb) where one of the leptons from the Z ££ decay is 
not reconstructed, result in the same final state signature 
as the Higgs boson in this search. Therefore, we search 
for these processes to validate our analysis methodology. 
The only change in the analysis is in the training of the 
final discriminant in events with two or three jets with 
one tight fo-tag or two 6-tags. We train using the WZ and 
ZZ diboson processes as signal while leaving the WW 
process as a background. The output of this discrimi- 
nant is used to measure the combined WZ and ZZ cross 
section by performing a maximum likelihood fit to data 
using signal plus background models, with maximization 
over the systematic uncertainties as described in detail in 
Sec. XII. The expected significance of the measurement 
using the MVA output is 1.8 s.d. We measure a cross 
section of 0.50 ± 0.34 (stat.) ± 0.36 (syst.) times the ex- 
pected SM cross section of 4.4 ± 0.3 pb. Figure 20 shows 
the MVA discriminant output for the diboson cross sec- 
tion {WZ + ZZ) with background-subtracted data and 
signal scaled to the best fit value. 
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FIG. 20: (color online) Final MVA discriminant output shown 
for the expected diboson signal and background-subtracted 
data rebinned as a function of log(5'/-B), after the maximum 
likelihood fit, summed over 6-tag channels. The error bars 
on data points represent the statistical uncertainty only. The 
post-fit systematic uncertainties are represented by the solid 
lines. The signal expectation is shown scaled to the best fit 
value. The inset gives an expanded view of the high log(5'/i3) 
region. 
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XII. UPPER LIMITS ON THE HIGGS BOSON 
PRODUCTION CROSS SECTION 

We derive upper limits on the Higgs boson production 
cross section multiplied by the corresponding branching 
fraction in units of the SM prediction. The limits are cal- 
culated using the modified frequentist CLs approach [75- 
77] , and the procedure is repeated for each assumed value 
of Mh- 

Two hypotheses arc considered: the background-only 
hypothesis (B), in which only background contributions 
are present, and the signal-plus-background (S+B) hy- 
pothesis in which both signal and background contribu- 
tions arc present. 

The limits are determined using the MVA output dis- 
tributions, together with their associated uncertainties, 
as inputs to the limit setting procedure. To preserve the 
stability of the limit derivation procedure in regions of 
small background statistics in the one tight 6-tag and all 
two 6-tags categories, the width of the bin at the largest 
MVA output value is adjusted by comparing the total 
background and signal-t-background expectations until 
the statistical significances for B and S+B are, respec- 
tively, greater than 3.6 and 5.0 s.d. from zero. The 
remaining part of the distribution is then divided into 
equally sized bins. In the zero 6-tags and one loose &-tag 
categories, the width of the bin at largest MVA output 
is set such that the relative statistical uncertainty on the 
signals plus background entries is less than 0.15. The 
remaining bins are distributed uniformly. The rebinning 
procedure is checked for potential biases in the determi- 
nation of the final limits, and no such bias is observed. 

We evaluate the compatibility of the data with 
the background-only and signal-|-background hypothe- 
ses. This is done using the log likelihood ratio (LLR), 
which is twice the negative logarithm of the ratio of 
the Poisson likelihoods, L, of the signal-|-background 
hypothesis to the background only hypothesis, LLR = 
~2\xy{Ls+b/Lb). 

Systematic uncertainties are included through nui- 
sance parameters that are assigned Gaussian probability 
distributions (priors). The signal and background pre- 
dictions are functions of the nuisance parameters. Each 
common source of systematic uncertainty (such as the 
uncertainties on predicted SM cross sections, identifica- 
tion efficiencies, and energy calibration, as described in 
Sec. X) is taken to be correlated across all channels ex- 
cept as otherwise noted in Sec. X. 

The inclusion of systematic uncertainties in the gen- 
eration of pscudocxpcriments has the effect of broaden- 
ing the expected LLR distributions and, thus, reducing 
the ability to resolve signal-like excesses. This degrada- 
tion can be partially reduced by performing a maximum 
likelihood fit to each pseudoexperiment (and data), once 
for the B hypothesis and once for the S-l-B hypothesis. 
The maximization is performed over the systematic un- 
certainties. The LLR is evaluated for each outcome us- 
ing the ratio of maximum likelihoods for the fit to each 



hypothesis. The resulting degradation of the limits due 
to systematic uncertainties is ~ 30% for searches in the 
vicinity of Mr = 125 GeV. 

The medians of the obtained LLR distributions for the 
B and S+B hypotheses for each tested mass are presented 
in Fig. 21. The corresponding ±1 s.d. and ±2 s.d. values 
for the background-only hypothesis at each mass point 
are represented by the shaded regions in the figure. The 
LLR values obtained from the data are also presented in 
the figure. 

The MVA discriminant distributions, for the Higgs 
boson mass point Mh = 125 GeV, after subtracting 
the total posterior background expectation are shown in 
Fig. 22. The signal expectation is shown scaled to the ob- 
served upper limit (described later) and the uncertainties 
in the background after the constrained fit are shown by 
the solid lines. 

Upper limits are calculated at 23 discrete values of 
the Higgs boson mass, spanning the range 90-200 GeV 
and spaced in increments of 5 GeV, by scaling the ex- 
pected signal contribution to the value at which it can 
be excluded at the 95% C.L. The expected limits are 
calculated from the background-only LLR distribution 
whereas the observed limits are quoted with respect to 
the LLR values measured in data. The expected and ob- 
served 95% C.L. upper limits results for the Higgs boson 
production cross section multiplied by the decay branch- 
ing fraction are shown, as a function of the Higgs boson 
mass Mh, in units of the SM prediction in Fig. 23. The 
values obtained for the expected and observed limit to 
SM ratios at each mass point arc listed in Table IV for 
all one-tight, two-loose, two-medium, and two-tight 6-tag 
subchannels together, for the two-jet and three-jet, zero 
and one loose 6-tag subchannels (all 6-tag categories for 
Mh > 150 GeV) together, the > 4-jet subchannels, and 
the combination of all subchannels. 



XIII. INTERPRETATIONS IN FOURTH 
GENERATION AND FERMIOPHOBIC HIGGS 
MODELS 

Extensions of the minimal electroweak symmetry 
breaking mechanism of the SM may be allowed, includ- 
ing models with a fourth generation of fermions or with 
a Higgs boson that has modified couplings to fermions, 
as in fermiophobic Higgs models (FHM). We interpret 
our results in these scenarios using the subchannels that 
are sensitive to 77 — )- WW decays: events with two or 
more jets and zero or one loose 6-tag for Mh < 150 GeV, 
extended to include pretag two- and three-jet events for 
Mh > 155 GeV. These are the first results for these 
models in the ^j/+jets final state. 

Previous results from the Tevatron Collider experi- 
ments in the context of a fourth generation of fermions 
set a limit on the Mh of 131 < Mh < 207 GeV [78]. 
The ATLAS [79] and CMS [80] collaborations exclude 
140 < Mh < 185 GeV and 144 < Mh < 207 GeV, re- 



24 



V(^ lv)+2, 3 jets with 1 tight+2 b-tags V(^ lv)+2, 3 jets with 0+1 loose b-tags 




100 120 140 160 180 200 100 120 140 160 180 200 

Mh (GeV) Mh (GeV) 



FIG. 21: (color online) The expected and observed log likelihood ratios as functions of the hypothesized Higgs boson mass Mh 
for the (a) electron and muon, two and three jets, one tight and two fe-tag channels; (b) electron and muon, two and three jets, 
zero and one loose 6-tag channels; (c) electron and muon, four or more jets, zero and one loose b-tag channels; (d) combination 
of all channels. The dashed red and black lines correspond to the median LLR of the signal+background and background-only 
hypotheses, respectively. The solid line corresponds to the LLR obtained from the data, and the shaded regions are the ±1 s.d. 
and ±2 s.d. values for the background-only hypothesis. 



TABLE IV: The expected and observed 95% C.L. limits, as a function of the Higgs boson mass Mh, presented as ratios of 
production cross section times branching fraction to the expected SM prediction. 
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spectively. Previous searches for the fermiophobic Higgs 
boson in — > 77 and H — > VV channels, with two 
leptons in the final state, were carried out at the LEP 
e+e- Colhder [81-84], by the CDF [85] and DO [86] 
Collaborations, and by the ATLAS [87] and CMS [88] 
Collaborations, with the most stringent limits being set 



by the CMS experiment where the excluded range is 
110 < Mh < 194 GeV. 

The Hgg coupling is enhanced in fourth-generation 
models, which leads to a higher rate of gg ^ H pro- 
duction and a larger decay width oi H gg than 
in the SM [89-92]. However, since H gg is loop- 
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FIG. 22: (color online) The MVA discriminant output dis- 
tribution minus the total background expectation for Mh ~ 
125 GeV rebinned as a function of log{S/B). The post-fit 
uncertainties are represented by the solid lines. The signal 
expectation is shown scaled to the best fit value. The inset 
gives an expanded view of the high log(S'/B) region. 



mediated, the H — > WW* decay mode dominates for 
Mh > 135 GeV, as in the SM. We consider two scenarios 
for the presence of a fourth generation. In the "low-mass" 
scenario, we assume a fourth-generation neutrino mass 
of 771ij4 ~ 80 GeV and a value for the fourth-generation 
charged Icpton mass oi rriii = 100 GeV, while in the 
"high-mass" scenario, we assume values for the fourth- 
generation neutrino and lepton masses of mi,^ = m£4 = 
1 TeV. Both scenarios set the fourth-generation quark 
masses to the values in Ref. [92]. After applying our se- 
lection criteria, the total expected signal for gg — > H pro- 
duction in the low-mass (high-mass) fourth-generation 
model is enhanced by a factor of 7.2 (7.5) over the SM 
production rate for Mh = 125 GeV. We only consider 
gluon fusion Higgs boson production, and wc set limits on 
a{gg — > i?) x B{H — > lyiy*). These Umits are compared 
with the predicted gg H production cross section re- 
sults from HDECAY [51], as shown in Fig. 24. We exclude 
the "low-mass" scenario for 150 < Mh < 188 GeV, and 
the "high-mass" scenario for 150 < Mh < 190 GcV. 

In the FHM, the Higgs boson does not couple to 
fermions at tree level but is otherwise SM-like. This sup- 
presses production via gluon fusion to a negligible rate 
and forbids direct decay to fermions. Production in asso- 
ciation with a vector boson or via vector boson fusion is 
allowed. For this interpretation, wc set the contribution 
from gg ^ H production to zero and scale the contri- 
butions from other production and decay mechanisms to 
reflect the predicted rate in the FHM. After applying our 
selection criteria, the total expected signal for vector bo- 
son fusion and VH — )- VWW production in the FHM is 
enhanced by a factor of 4.2 over the SM production rate 
for Mh = 125 GeV. The expected and observed cross 
section times branching fraction limits are compared to 
the FHM predictions in Fig. 25. 



XIV. SUMMARY 

We have presented a search for SM Higgs boson pro- 
duction in lepton + + jets final states with a dataset 
corresponding to 9.7 fb~^ of integrated luminosity col- 
lected with the DO detector. The search is sensitive 
to VH Vbb, H ^ WW* ii^jj, and WH 
WWW* — >■ ivjjjj production and decay, and super- 
sedes previous VH — J> Vbb and H WW* Ivjj 
searches published by DO. To maximize our signal sen- 
sitivity, we subdivide the dataset into 36 independent 
subchannels according to lepton flavor, jet multiplicity, 
and the number and quality of 6-taggcd jets and apply 
multivariate analysis techniques to further discriminate 
between signal and background. We test our method by 
examining SM WZ and ZZ production with Z ^ bb de- 
cay and find production rates consistent with the SM 
prediction. We observe no significant excess over the 
background prediction as expected from the amplitude 
of a 125 GeV SM Higgs boson signal, given the sensi- 
tivity of this single channel. Significance is achieved by 
combining this channel with the other low mass chan- 
nels analyzed at the Tevatron [13], while here we set 
95% C.L. upper limits on the Higgs boson production 
cross section for masses between 90 and 200 GeV. For 
Mh = 125 GeV, the observed (expected) upper limit 
is 5.8 (4.7) times the SM prediction. We interpret the 
data also in models with fourth generation fermions, or 
a fcrmiophobic Higgs boson. In these interpretations, we 
exclude 150 < Mh < 188(190) GeV in the "low-mass" 
("high-mass") fourth generation fermion scenario, and 
provide 95% C.L limits on the production cross section 
in the fcrmiophobic model. 
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FIG. 23: (color online) The expected and observed 95% C.L. upper limits on SM Higgs boson production for the (a) electron 
and muon, two and three jets, one tight and two &-tag channels; (b) electron and muon, two and three jets, zero and one loose 
6-tag channels {Mr < 150 GeV) and pretag channels (A/h > 155 GeV); (c) electron and muon, four or more jets, zero and one 
loose 6-tag channels; (d) combination of all channels. The limits are presented as ratios to the expected SM prediction. The 
dashed line corresponds to the expected limit, and the solid line corresponds to the limit observed in data. The shaded regions 
are the ±1 s.d. and ±2 s.d. values for the expected limit. 
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FIG. 24: (color online) The expected and observed 95% C.L. 
upper limits on (j{gg — H) x B{H — >■ VKVK) compared to the 
prediction from the fourth-generation fermion model. 



FIG. 25: (color online) The expected and observed 95% C.L. 
upper limits on fermiophobic Higgs boson production. 
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Appendix A: Multivariate Discriminator Input 
Variables 

The multivariate discriminators used in this search use 
input variables from five general categories: final state 
particle information, as measured in the DO detector; 
kinematics of reconstructed objects, such as W boson 
candidates reconstructed from the leptonic or hadronic 
decay products; angular distributions between final state 
particles and reconstructed objects; topological variables 
that examine the net properties of all final state parti- 
cles in an event; and special variables focused on dis- 
criminating Higgs boson candidate events from specific 
backgrounds. Certain multivariate discriminants trained 
to separate a Higgs boson signal from a specific back- 
ground are also used as inputs for a final discriminant 
that is trained to separate the Higgs boson signal from 
all backgrounds. 

Individual input variables are described in detail be- 
low. In the descriptions, i refers to the electron or muon 
in a selected event, u refers to the neutrino candidate, 
and jn refers to jets as ordered by pr where ji is the 
jet with highest pt- The pz of the neutrino candidate 
is estimated by constraining the charged lepton and the 
neutrino system to the mass of the W boson and choosing 
the lowest magnitude solution. 

Input variable lists for each multivariate discriminant 
appear in Tables V-XVIII. The ranking by importance of 
the variables is determined in the BDT by counting how 
often the variables are used to split decision tree nodes, 
and by weighting each split occurrence by the separation 
gain-squared it has achieved and by the number of events 
in the node [33]. In the RF, the importance of variables 
are estimated after training in an independent sample of 
validation events. These events are run through the RF, 
once for each variable used. On each pass the class of 
each event is randomized whenever the variable under 
test is encountered and the change in the quadratic loss 
figure of merit is estimated: 

^events 

wgUieventf"'' ~ RF{evenU) f 
FOM = (Al) 

i^events 

^ wgU 

1=1 

where wgti is an event weight, eventf°''^'^=\ for signal, 
for background, and RF[eventi) is the output of the RF 
classifier for a given event. Whenever the RF makes an 
incorrect assignment for an event the FOM increases in 
value. In this test the assignments are randomized for 
one variable at a time, effectively removing the predic- 
tive power of that variable, and the FOM will increase 
more when more powerful variables are removed in this 
manner. 

The input variable distributions arc defined as follows: 



1. Final State Particle Information 

• E{ji): Energy of the leading jet 

• p-l^ : Pt of the leading jet 

• pi^: Pt of the second leading jet 

• Pt ■ Pt of the third leading jet 

• X rj^: Product of the lepton charge and its pscu- 
dorapidity 

• : Smaller absolute value solution for pz of the 
reconstructed neutrino, reconstructed with the as- 
sumption all ^T is originating from W boson 

• p'^^ : Larger absolute value solution for pz of the 
reconstructed neutrino 

• ;^T- Missing pT as determined from charged parti- 
cle tracks in central tracking detector 

• ^T^^: Scaled ^t is defined as 

• _^^'^: ^T significance, a measure of the consistency 
of the observed ^t with respect to zero ^t, ac- 
counting for the uncertainty in the calorimeter ob- 
jects that contribute to ^t 

• ?^T^'^: significance, a measure of the consistency 
of the observed ^t with respect to zero ^t, ac- 
counting for the uncertainty in the charged particle 
tracks that contribute to ^t 

9 q^x rj^^ : Product of the the lepton charge and pseu- 
dorapidity of the second leading jet 

• q^ X rj^^ : Product of the the lepton charge and pseu- 
dorapidity of the third leading jet 

• : Averaged 6-jet identification output for the 
highest energy 6-tagged jets 

2. Kinematics of Reconstructed Objects 

• ?7iji2 : Invariant mass of the leading and second 
leading jets 

• mPr^'^ : Transverse mass of the leading and second 
leading jets 

• mjj23 : Invariant mass of the leading, second lead- 
ing, and third leading jets 

• ?nji23_j : Invariant mass of the leading, second lead- 
ing, third leading, and fourth leading jets 

• Apt{1,^t)- Scalar difference: \p^ ~ ^t\ 

• ApT{ji,j2)- scalar difference, ^ Pt 
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J2PTijijj2,£)'- Scalar sum of the pr of the two 
leading jets and the lepton 

Apt{£,^t)/Pt^ ■ Ratio of the scalar difference be- 
tween pfp and the , to p^ 

niax{p^ , ^t) / Pt' ■ Ratio of the niax(p^,^T) to p}p 
nim{p^ , ^t) / Pt^ ■ Ratio of the niin(p^,^T) to p^ 
pf/{pf^ +^t)- Ratio of the pf to + 
ApTiWJ): \p^~p't\ 
Apt{W,^t): \p^~^t\ 

p-!^^^: pt of the system consisting of the leading, 
second leading, and third leading jets 

ApT{j2T 325)'- scalar Apx between the second lead- 
ing jet and the system consisting of the second lead- 
ing and third leading jets 

2 

^ p!^ : Scalar sum of the pT of the two leading jets, 
1=1 

Pt +p't 

4 

^ : Scalar sum of the pt of the leading, second 

1=1 

leading, third leading and fourth leading jets 
2 

Pt^ I Til Pt'- Ratio of the pt of the leading and 

4=1 

second leading jet system to the scalar sum of the 
Pt of the two leading jets, 

3 

Pt^ I TjPt'- Ratio of the pt of the system consist- 

i=2 

ing of the second leading and third leading jets to 
the scalar sum of the pT of the second leading and 
third leading jets 

Recoil(p^^^): Recoil pt of the first and second lead- 
ing jet system 

rrij-^^t'- Invariant mass of the dijet system and the 
lepton 

M^: Transverse mass of the iiy system 

p^ : pt of the iv system 

pip + ]^T'- Scalar sum of pip and ]^t 

Recoil(p^): pt of the tv system with respect to 
the thrust vector, (. — v 

miiijx '■ Invariant mass of the system consisting of 
the charged lepton, reconstructed neutrino (assum- 
ing p'2^), and leading jet 

nii^j^ : Invariant mass of the system consisting of 
the charged lepton, reconstructed neutrino (assum- 
ing p'^), and second leading jet 



• p^^"'^: Pt of the system consisting of the charged 
lepton, reconstructed neutrino (assuming ), and 
second leading jet 

• Pt+ ^t + Pt ■ Scalar sum of the pT of the charged 
lepton, ^T, and second leading jet 

• 'me,uji2 ■ Invariant mass of the charged lepton, recon- 
structed neutrino (assuming ) j and two leading 
jets 

• rrip''^^: Transverse mass of the charged lepton, ^t, 
and two leading jets 

• me^j-^2{pz{v) = 0): Invariant mass of the charged 
lepton, reconstructed neutrino (assuming ), and 
two leading jets, with the assumption that p'^ = 

2 

• Y.PT{t,^T, 31,32)- Scalar sum, + I] Pt 

1=1 

• I]pt(^,^t, ji, J2, J3,J4): Scalar sum, pip + $t + 
tp^T 

i=l 

• H: Helicity is defined for an object A, coming 
from the decay of object C via C -> AB, as 
arccos0^c[(^-^)/(|C| x \A\)] 

• "^(^12, ji): Helicity of the leading jet in the dijet 
system, calculated in the laboratory frame 

• V: Velocity is defined for an object C AB as 
- ln{l - {1 - 4 X [{ml + ml)/ml]'/^rf'} 

• V(ji2): Velocity of the dijet system 

• Vji3- Velocity of the system consisting of the lead- 
ing and third leading jets 

• T: Twist is arctan(A0/Ary) 

• Tijri)' Twist of the dijet system 

• 7523' Twist of the system consisting of the second 
leading and third leading jets 

• Tw~^iu- Twist of the £iy system 

• W: Width of a jet in {-q, (f) space defined as 
^J^ni+^u-, where -q^ and (f) are the pt weighted 
RMS T] and (j) of energy deposits around the jet 
centroid. 

• yVj3'- Width of the third leading jet 
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3. Angular Distributions 

• Ar]{ji , 72): Separation in 77 between the two leading 
jets, \r]j, -TjjJ 

• max I A7y(ji2, {ji or j2})h Maximum Ar) between 
the dijct system and the leading or second leading 
jet 

• , j2 ) : Separation in <j) between the two leading 
jets, - <j>jj 

• A72.(ji,j2): Angular separation in (rj^ip) space be- 
tween the two leading jets 

• min A7?,(ji2, {ii or J2}): Minimum angular separa- 
tion in (77, 0) space between the dijet system and 
the leading or second leading jet 

• A(/)(ji2, ji): - where ^j^^ is the of the 
dijet system 

• A(f)(ji, js): Separation in between the first and 
third leading jets, \(j>j-^ — (fij^ \ 

• A7^(j3, jia): ATZ between the third leading jet and 
the system consisting of the leading and third lead- 
ing jets 

• A(/)(73, J23): A(f) between the third leading jet and 
the system consisting of the second leading and 
third leading jets 

• j2): 3D angle between the two leading jets 

• cos J2)cm: Cosine of the angle between the 
two leading jets in the center of mass (CM) of the 
dijet system 

• Z[^,bis(ji, ^2)]: 3D angle between the charged lep- 
ton and the bisector of the dijet system 

• A0[VF,bis(ji, J2)]: Signed Ai^ between the iiy sys- 
tem and the bisector of the dijet system 

• |A77(£,ji)|: Separation in 77 between the lepton and 
the leading jet, [77^ — \ 

• |A77(^,j2)h Separation in 77 between the lepton and 
the second leading jet, [77^ — 77-'^ | 

• |A77(£,j3)|: Separation in 77 between the lepton and 
the third leading jet, \r]^ — r]^^\ 

• max |A77(^, {ji or .72})h Maximum A77 between the 
charged lepton and the leading or second leading 
jet 

• ATZ{£, ji): ATZ between the charged lepton and the 
leading jet 

• ATZ{£, J2): ATZ between the charged lepton and the 
second leading jet 



• A72.(€, js): ATZ between the charged lepton and the 
third leading jet 

• A(/)(^T, ji): A0 between the and the leading 
jet 

• ATZ{iy, ji): ATZ between the reconstructed neutrino 
(assuming p^^), and the leading jet 

• i[mn[ATZ{i>, {ji or J2})]: Minimum ATZ between 
the reconstructed neutrino (assuming p'^ ) , and the 
leading or second leading jet 

• A7?,(€, light jet): ATZ between the charged lepton 
and the leading non-&-tagged jet 

• Z(€, J12): 3D angle between the charged lepton and 
the dijet system 

• A77(£, I'): Separation in rj between the lepton and 
the reconstructed neutrino (assuming ), \'rf^i]'^\ 

• max |A77(VF, {i or Maximum A77 between the 

system and charged lepton or reconstructed neu- 
trino (assuming ) 

• A(f){l,]^T)'- 4> angle between the lepton and ]^t- 

• max |A(^(W^, {£ or ^T})h Maximum A0 between 
the iv system and the charged lepton or 

• min I A0(iy, {£ or ^T})h Minimum A0 between 
the system and the charged lepton or 

• ATZ{i^ v): ATZ between the charged lepton and the 
reconstructed neutrino (assuming ) 

• max[A7?,(VF, or I/})]: Maximum ATZ between 
the system and the charged lepton or recon- 
structed neutrino (assuming p^^ ) 

• mva[ATZ{W, {£ or ly})]: Minimum ATZ between the 

system and the charged lepton or reconstructed 
neutrino (assuming ) 

• |A7;(VF, £)|: A77 between the system and the 
charged lepton 

• ATZ{W,£): ATZ between the system and the 
charged lepton 

• ATZ{W, v): ATZ between the £v system and the re- 
constructed neutrino (assuming ) 

• /L{(.^v): 3D angle between the charged lepton and 
the reconstructed neutrino (assuming ) 

• cos9{t)i^cyi'- Cosine of the angle between the 
charged lepton and the proton beam axis in the 
CM of system 

• cos 9{l): Cosine of the angle between the charged 
lepton and the proton beam axis in the detector 
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• |A77(Ty, j2)|: Ary between system and the second 
leading jet 

• A(j){W,j2)- between £v system and the second 
leading jet 

• ATZ{W,j2)- AT?, between system and the second 
leading jet 

• |A?7(M^, ji2)|: A77 between £1/ system and the dijct 
system 

• ^{ji,j2)}iCM- 3D angle between the two leading 
jets in the H WW Ivjj CM frame (HCM) 

• cos[Z(ji,.^zy)]HCM: Cosine of the 3D angle between 
the leading jet and Iv system in the HCM 

• cos[Z(ji, £)]hcm: Cosine of the 3D angle between 
the charged lepton and the leading jet in the HCM 

• cos[Z(jf (W^ ^ £i/)HCM)]. Cosine of the 3D 
angle between the leading jet in energy in the CM of 
the dijet system and £v system in the H — > WW — > 
ivjj CM frame; jet energy is calculated in the H — >■ 
WW tujj CM frame 

• cos[/L{£i^cMt{W — > £v)yicm)\'- Cosine of the 3D 
angle between the charged lepton in the system 
CM and Iv system in the H WW -> ivjj CM 
frame 

• cos.m£iuCM-M,{W iv)i^jjCM-4])]- Cosine of 
the 3D angle between the charged lepton in the £v 
system CM and £v system in the H — > WW — > ivjj 
CM frame for V{-^ jj)H{-^ WW — >■ ivjj) candi- 
date events; jet energy is calculated in the H ^• 
WW £vjj CM frame 

• cos(Z(ji,£i^))HCM;4j: Cosine of the 3D angle be- 
tween the leading jet and iu system in the H — >■ 
WW — >■ ivjj CM frame frame for y(— > jj)H{-^ 
WW — > ivjj) candidate events 

• cos{9*): 9* — Z(W^, incoming w- type quark) in 
HCM frame [93] 

• cos(x*): X* — ^(^,spinn/) in iv system CM 
frame [93] 

4. Topological Variables 

• A: Aplanarity is 3A3/2 where A3 is the small- 
est eigenvalue of the normalized momentum tensor 

= (E.P?pf)/(EJpin > where a,/3 = 1,2,3 
correspond to the x,y,z momentum components, 
and i runs over selected objects. Without argu- 
ments, it is calculated for all visible objects 

• A{ijij2)'- A calculated for the charged lepton, and 
leading and second leading jets 



• Aivi): A calculated for the charged lepton, recon- 
structed neutrino (assuming p^^), and all selected 
jets 

• C: Centrality is (^^Pt) l(^i\'Pi\)j where i runs 
over i and all jets 

• S: Sphericity is 3(A2 -I- A3)/2 where A3 (A2) is 
the smallest (second-smallest) eigenvalue of the 
normalized momentum tensor described under A. 
Without arguments, it is calculated for all visible 
objects 

• S{ijij2)- S calculated for the charged lepton, and 
leading and second leading jets 

• S{iv2jij2)' S calculated for the charged lepton, 
reconstructed neutrino (assuming p^^), and leading 
and second leading jets 

• S{vi): S calculated for the charged lepton, recon- 
structed neutrino (assuming p^^), and all selected 
jets 

• p^^^ : Magnitude of the vector sum of the px of the 
visible particles 

• J^Ipt)^^^ ■ Scalar sum of the pt of the visible par- 
ticles 

• Kf"": A7^(jl, J2) X Ei^/{]^T + -^t), where Et is 
the transverse energy 

• bis(ji2,!^): Scalar product of the bisector of the 
dijet system and the vector, i.e. bis(ji, ^2) • 

• m**'^'": Mass asymmetry between iv system and 
the dijet system: (m^^ — ^ji2) / {"^^v + ™ii2) 

• {p't ~Pt^ ) I {Pt +Pt^)' Pt asymmetry between iv 
system and the dijet system 

• Pt IPt''- Ratio of p^ to p^'- 

2 

• iPT+^T-'T,PT)/T,PTii,^T,ji,j2): X^PTasym- 

1=1 

metry between iv system and the dijet system 

2 2 

• {p^j. +^t)/EPt- R-atio of +^TtoJ2 Pt 

1=1 1=1 

• pp'yEPT{i,^T,ji,j2y- Ratio of the pp'^ to 

Y^PT{i,^T,jl,j2) 

• SIG{ji2, ji)- Based on the pull variables de- 
scribed in Ref. [94]. Sigma, SIQ, of the dijet 
system with respect to the leading jet defined as 

pi^ATZin,j2)/j:PT 

8=1 
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• max[5T^(ji2, {ji or ^2})]: Maximum STQ of 
the leading or second leading jet defined as 

2 

^'T''''0'l'J2)A7^(jl,J2)/ Y.Pt '^ith respect to the 
1=1 

dijet system 

• min[5Xf/(j'i2, {j'l or J2})]: MinimmTi SIQ of 

the leading or second leading jet defined as 
2 

PT""(ji'J2)A7e(ji,j2)/ E Pt with respect to the 
dijet system 

• SIQ{W,t): SIQ of the £v system with respect to 
the lepton, defined as (A7^(€, v) x ^t)I(j>t + V't) 

• max[5X5(W, or i^})]: Maximum SXQ of the lep- 
ton or defined as max(p^,^T) x ATZ{i, v)/{p^ + 
^t) with respect to the system 

• T[m\[SIQ{W,{i or v})]: Minimum SIQ of the lep- 
ton or defined as min(p^, ^t) x A7?.(^, v)/{p^rp 4- 
]^t) with respect to the system 

• STQ-^cts{A) : Sigm defined for an object A, with 
respect to all jets in an event, as 

Ejotsbf xA7^(AJet)]/Ejet.b^^') 

• 5It/jets(ji23): ^Itjjcts of the system consisting of 
the leading, second leading and third leading jets 

• iSIC/jcts(ji234): SIQjets of the system consisting 
of the leading, second leading, third leading and 
fourth leading jets 

• SIQjcts{i)- SIQ jets of the charged lepton 

• SIM: Similarity is defined for two objects, A and 
B, as min(pT)^ x A7?.^/(X]pt)^, where min(pT) is 
the minimum px of the objects A and B, A72. is the 
angular separation in (r;, 0) space between objects 
A and B, and J2pt is the scalar sum of the prs of 
A and B 

• SIM {£, v) : Similarity of the charged lepton and 
the reconstructed neutrino (assuming ) 



• MVAMj(T^-ff): Output of the multivariate discrim- 
inant trained to distinguish WH Ivhh from the 
MJ background 



5. Discriminants to Separate Higgs Boson Events 
from a Specific Background 

• MVAtt: Output of the multivariate discriminant 
trained against tt and single top quark backgrounds 

• MVAyj: Output of the multivariate discriminant 
trained against V^-|-jets backgrounds 

• MVAvv- Output of the multivariate discriminant 
trained against diboson backgrounds 

• MVAMj(i^ YV): Output of the multivariate 
discriminant trained to distinguish B — > WW — > 
ivjj from the MJ background 
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TABLE V: Input variables for the MVAm3{VH) and 
MVAm.j(J^ — >■ VV) discriminants. Two discriminants were 
trained to reject MJ events, one trained using VH — ^ ivbb 
events as a signal and the second using H ^ VV ^ liyjj 
events as signal. Both discriminants use the same list of in- 
puts. Variables are listed by their importance in the MVA. 



MVAmj Input Variables 

Pt 

V(jl2) 
^Asym 

c 

„vis 
Pt 

max|A??(^, {ji or j2})| 



TABLE VI: Table of input variables for the final signal dis- 
criminant for the WH — >■ Ivbh channel. Variables are listed 
by their rank of importance when used in the two tight b- 
tagged (2T), two medium 6-tagged (2M), two loose fe-tagged 
(2L), and one tight 6-tagged (IT) categories. 



Variable 


2T 


2M 


2L 


IT 


MVAmj (Vi/) 


1 


1 






mbib2 


2 


4 


3 


1 


p¥/{PT + m 


3 


6 


4 


2 




4 


13 


1 


4 


cos(x*) 


5 


3 






max|A77(f, {ji or j2})\ 


6 


11 


2 


3 


q X n 


7 


2 


6 


6 


A7^(^,il) 


8 


5 






min[52:C/(jii2, {j'l or ^2})] 


9 


15 


9 


5 




10 


7 


11 


9 


V(il2) 


11 


12 


7 


11 


cos(e*) 


12 


10 








13 


16 


12 


13 




14 


14 






C 


15 


8 


8 


10 


uprr' 


16 


9 






^Asym 






5 


8 


A 






10 


12 


Pt 






13 


7 
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TABLE VII: Table of input variables for the final signal dis- 
criminant for the H — >■ WW — >■ evjj channel for two jets 
events. Variables are by their rank of importance when used 
in the zero fo-tags (OT) and one loose 6-tag (IT) categories for 
Mh < 150 GeV. 





U 1 


1 ' 1 
1 1 


2 

(p^+^t)/Ep't 

i—l 


1 


1 


\Ari{W,£)\ 


2 






3 


4 




4 


6 


\Arj{W,j2)\ 


5 




max A77(ji2, {ji or j2}) 


6 


10 


A{£jij2) 


7 




A7^(^,il) 


8 






9 


15 


H(ll2, 71 ) 


10 




max|A0(VF,{^ or ^t})\ 


11 






12 




^Pt {31,32) 


13 






14 


7 


|A77(^,ji)| 




2 






3 


|A77(^,J2)| 




5 


V(il2) 




8 






9 






11 






12 


Apt{W,Pt) 




13 


Pt+Pt 




14 


Y^PT{jl,j2,£) 




16 


tp^ 




17 



TABLE VIII: Table of input variables for the final signal dis- 
criminant for the H — >■ WW — >■ evjj channel for three jets 
events. Variables are listed by their rank of importance when 
used in the zero fe-tags (OT) and one loose fe-tag (IT) cate- 
gories for Mh < 150 GeV. 



Variable 


U i 


i i 


A/TA/'A (Xl U\ 

MVAmjI.I' -n j 


1 
i 


n 
D 


2 




1 K 

io 


{p't+Pt)/T.P^ 


3 


2 


An{£,ji) 


4 




\^V{W,j2)\ 


5 




\Arj(W,£)\ 


6 




max 1 Ay7(7i2, In or 72|)| 


7 


8 


A{£jij2) 


8 




E(pt)'''' 


9 


16 


max|A(/)(M^,{^ or ^t})1 


10 




Apt ( 71 ^ 72 ) 


11 




'T-VJ12 1 J17 


12 




A<I>{£,^t) 


13 




Pt 


14 




MWAm}{H VV) 




1 


|A77(^,i2)l 




3 


|A?7(^,ii)| 




4 


COs6l(^)toCM 




5 






7 


V(il2) 




9 






10 






11 


Apt{W,^t) 




12 


Pt+^t 




13 


Y^PT{jl,j2,£) 




14 


tp^ 




17 



1=1 



36 



TABLE IX: Table of input variables for the final signal dis- 
criminant for the H — >■ WW —> ^ivjj channel for two jets 
events. Variables are listed by their rank of importance when 
used in the the zero 6-tags (OT) and one loose 6-tag (IT) 
categories for Mh < 150 GeV. 





OT 


IT 


A,^[W^,bis(ji,j2)] 


1 




\AviW,l)\ 


2 




'^312 


3 


7 




4 


16 


min[A7^(W', {£ or i^})] 


5 




A7^(jl,j2) 


6 


14 


SIM{£,u) 


7 




2 

1^ Pt 

1=1 





A 
4 


min(pr,_^T)/p?' 


9 




ApTii,^T) 


10 




A7^(//,Jl) 


11 




A<^(^T,Jl) 


12 




sig{ji2,ji) 


13 






14 




ApT [31,32) 


15 




p^ 


16 


6 


\^v{£,ji)\ 




1 


W0i2,ii) 




2 


\^vie,j2)\ 




3 


COs6'(^)toCM 




5 


A7^(W,j2) 




6 


|Ar,(W,ii2)| 




8 






9 


bis(ji2, 1^) 




10 


ApTiW,£) 




11 






12 


Tw^iv 




13 


A</>(ji,i2) 




15 






17 


J2pt{3i,32,£) 




18 


ATZ{£, u) 




19 


A<P{£,Pt) 




20 


p¥ 




21 



TABLE X: Table of input variables for the final signal discrim- 
inant for the H — >■ WW — >■ channel for three jets events. 
Variables are listed by their rank of importance when used in 
the zero fe-tags (OT) and one loose fe-tag (IT) categories for 
Mh < 150 GeV. 



Variable 


OT 


IT 




1 




An{£,j2) 


2 




\AviW,£)\ 


3 




\^vi£,3i)\ 


4 


2 


COs6'(ji, j2)cM 


5 


22 


A0(jl,j3) 


6 




V(il2) 


7 


21 


Recoil(p?') 


8 


16 


A 


9 






10 






11 


18 




12 




C0s[Z(jl,^!/)]HCM 


13 






14 




A'^(j3, J23) 




1 






3 


A<t>{W,j2) 




4 


JT SC 




5 


S{£hj2) 




6 


^3123 




7 


ATZ{v,n) 




8 


ATZ{£, light jet) 




9 


Ar?(ji, j2) 




10 


^Asym 




11 


Pt'/{Pt+^t) 




12 


SIQ jctsijl23) 




13 


C 




14 


sig{ji2,ji) 




15 


cos6'(^)«^CM 




17 






19 


mti,n^{pz{v) = 0) 




20 






23 






24 


An{£..v) 




25 
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TABLE XI: Table of input variables for the final signal dis- 
criminant for the H — >■ WW — >■ evjj channel for two jets 
events in the pretag category for Mr > 155 GeV. Variables 



TABLE XIL Table of input variables for the final signal dis- 
criminant for the H — > WW — ^ ei^jj channel for three jets 
events in the pretag category for Mr > 155 GeV. Variables 



■ importance in the MVA. 


are listed by their importance in the MVA. 


MVA Input Variables 


MVA Input Variables 






Recoil(p^) 


Recoil(p?') 


^312 


MVAmj(V//) 


j„Asym 


uia.^\Act>{W,{£ or ^t})\ 


sig(w,£) 


Sig{W,£) 


MYAmj{VH) 


^Asym 


\Av{W,£)\ 


\Av{W,£)\ 


E{h) 


^312 


An{£,j,) 


An{£,n) 


inaji\A(j>{W,{£ or ^t})\ 


E{ji) 




Ah, 32) 


ATi{W,£) 


C0s[Z(jl,/!/)]HCM 


Recoil(p^" ) 


W(ii2,ii) 


COs[Z(ji,^l/)]HCM 


Recoil(p^i" ) 


^0'l2,Jl) 


A7^(^, light jet) 


A7^(£, light jet) 
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TABLE XIII: Table of input variables for the final signal dis- 
criminant for the H — >■ WW — >■ fJ^i'jj channel for two jets 
events in the pretag category for Mr > 155 GeV. Variables 
are listed by their importance in the MVA. 

MVA Input Variables 
VO12) 

A<p[W,hmiji,j2)] 
min[A7^(^^, {ji or 72})] 

COs[Z(jl,^)]HCM 

A7^(il,i2) 

1=1 

min A7^(jl2, {ji or ji}) 

J2PT{£,^T,jl,j2) 

min|A(^(W,{^ or ^t})| 
min[A7^(W', {I or u})] 

A<^(ji2,ii) 

max[Sig{W, {I or i^})] 
AprOi, J2) 
max(p|,,^T)/pT 
5XC/(ji2, Ji) 

min(pr, ^t)/p?" 



TABLE XIV: Table of input variables for the final signal dis- 
criminant for the H — >■ WW — >■ ^J.vjj channel for three jets 
events in the pretag category for Mh > 155 GeV. Variables 
are listed by their importance in the MVA. 

MVA Input Variables 

A7^(,n,i2) 

{p'¥ -p^')/{pt +pr) 

^O'l, j2)HCM 
A</)[W,bis(ji,j2)] 

An{e,ji) 

mi^j^^ipzir^) = 0) 

cos d{i)l^CM 
COS^O'l, j2)cM 

minlSigiW, {£ or //})] 

p¥/iPT+m 

sig{ji2,ji) 

'5I5jets(il23) 

r(ji2) 

pi^ 
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TABLE XV: Table of input variables for the MVAtt discrim- 
inant for the H — >■ WW — >■ tvjjjj channel. Variables are 
listed by their importance in the MVA. 



TABLE XVI: Table of input variables for the MVAv.j dis- 
criminant for the H — >■ WW — >■ ii^jjjj channel. Variables are 
listed by their importance in the MVA. 



MVA Input Variables 
I]PT(^.-^T,ii, J2,i3, J4) 

A7^(Ai3) 

^(jl, j2)HCM 

{p¥ -p^')/{p¥ +pr-) 

E{1) 

^Asym 

PT''VEPT(^,^T,il,i2) 

S 

\^ri{W,l)\ 

COs(Z(jl'i,^i^))HCM;4i 

cos[Z(if=^^,(W^/i/)f^c^) 
A0[Vy,bis(ii,i2)] 

ApT (j2,j23) 
A7^(j3,jl3) 

ATLiW, t) 

Z(£, v) 
Pz- 

RecoiKp?") 
Pt/p^' 

min|A<^(iy,{£ or ^t})| 



MVA Input Variables 



W,3 

COSe(^)toCM 

A7^(^,j3) 

"^23 

max|A7?(£,{ji or ^2})! 

{p'T+]?T-j:P^)/Y.PT{i,pT,juj2) 

max|A77(W^,{^ or u})\ 
A.TL{£,v) 
A7?(ii,i2) 
A7^(VF,^) 

COs[Z(^£i,CA/;4j, (W — >■ ^!^)<'^jjCMi4j)] 
C 

A7^(Ai2) 

COs[Z(^toCA/, {W -)■ ^!/)hCm)] 
i = l 

|Ar?(W,£)| 

w 
Pt 

A{vi) 

m 

Recoil(p?^) 

P^/p't'' 
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TABLE XVII: Table of input variables for the MVAvv dis- 
criminant for the H — >■ WW — >■ tvjjjj channel. Variables are 
listed by their importance in the MVA. 



MVA Input Variables 



COs6l(^)toCM 



CM 



{W 



-iHCM-i 



Ep'T)/EPT(^,^T,jl,j2) 

1=1 

maxjA77(£, {ji or j2})| 
q X r)^ 
X r;^^ 

{p¥ ^pr)/ip¥ +Pr) 

^O'li j2)HCM 

Aniw,i) 

Pp'Vi:PT{l,^T,jl,j2) 

S1M{£, v) 

A7^(Ai2) 

coa[A{£i^CM-Aj,{W -)■ iy)ivjjCM-Aj)] 

5IC/jets(il23) 

5IC/jots(il234) 

Z[£,bis(jl,j2)] 

A7^(j3, jis) 

i=l 

max[A7^(VF, or u})] 



TABLE XVIII: Table of input variables for the final signal 
discriminant for the H — >■ WW — >■ ii^jjjj channel. Variables 
are listed by their rank of importance when used in the zero b- 
tags (OT) and one loose 6-tag (IL) categories for each charged 
lepton type. 





c U -L 






MVAtt 


1 


r 




3 


MVAwT('ff — > VV) 


2 


7 


4 9 


MVAvj 


3 


1 


1 1 


MVAvv 


4 


8 


7 2 




5 


12 


14 




6 














P'r/T.PT 

i = 2 


7 




9 5 


77.Sig 

Pt 






1 /I 
14 




M V Am.j [y ^ } 





Q 


LA 14 


2 


10 


6 


16 




11 






q X 77 


12 




3 8 


COs[Z(^f^CM, {W — >■ ^i^)HCM)] 


13 








14 








15 






A7^(ly,^) 


16 




13 




17 


11 




COS e{£) 


18 


9 




Pt'^ 




2 








4 




SXQje^s (7123) 




10 




g*^ X 77^1 




13 




|A77(W,^)| 




15 


13 


max A77(VK, or u})\ 




16 




cos 0(71, 72)cM 






2 12 


c 






6 


|Ar7(^,j3)| 






8 


COs[Z(ji,^!/)]HCM 






10 








11 








15 


Z[Abis(ii,i2)] 






16 


C0s6'(^)f^CM 






17 








18 


A77(£,/.) 






19 








20 



A7^(^,i3) 

COs[Z(ji,^)]HCM 

mi^j^^{pz{v) = 0) 
5(^jii2) 

max[A7^(VK, or i^})] 
max[5IC;(ji2, {ji or ,72})] 
Recoil (p?") 



4 
6 

7 
10 
11 
15 
17 



